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Abstract 

The integrals of motion of the classical two dimensional superinte- 
grable systems with quadratic integrals of motion close in a restrained 
quadratic Poisson algebra, whose the general form is investigated. 
Each classical superintegrable problem has a quantum counterpart, 
a quantum superintegrable system. The quadratic Poisson algebra 
is deformed to a quantum associative algebra, the finite dimensional 
representations of this algebra are calculated by using a deformed 
parafermion oscillator technique. It is shown that, the finite dimen- 
sional representations of the quadratic algebra are determined by the 
energy eigenvalues of the superintegrable system. The calculation of 
energy eigenvalues is reduced to the solution of algebraic equations, 
which are universal for all two dimensional superintegrable systems 
with quadratic integrals of motion. 
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I Introduction 



In classical mechanics, integrable system is a system possessing more con- 
stants of motion in addition to the energy. A comprehensive review of the 
two-dimensional integrable classical systems is given by Hietarinta [[J , where 
the space was assumed to be flat. The case of non flat space is under current 
investigation [|. An interesting subset of the totality of integrable systems 
is the set of systems, which possess a maximum number of integrals, these 
systems are termed as superintegrable ones. The Coulomb and the harmonic 
oscillator potentials are the most familiar classical superintegrable systems, 
whose their quantum counterpart has nice symmetry properties, which are 
described by the so(N + 1) and su(N) Lie algebras respectively. 

The Hamiltonian of a classical system is generally a quadratic function of 
the momenta. In the case of the flat space, all the known two dimensional su- 
perintegrable systems with quadratic integrals of motion are simultaneously 
separable in more than two orthogonal coordinate systems M . The integrals 
of motion of a two dimensional superintegrable system in flat space close in 
a restrained classical Poisson algebra [|, ||, [Tj. The study of the quadratic 
Poisson algebras is a matter under investigation, related to several branches 
of physics as: the solution of the classical Yang - Baxter equation ||, the 
two dimensional superintegrable systems in flat space or on the sphere f?|, 
the statistics or the case of "exactly solvable" classical problems [10]. 

The quantization of a classical integrable system corresponds generally 
to a quantum integrable system. The mechanism of quantum deformation 
of a classical system to a quantum one is not fully understood. Initially the 
problem of quantization of classical superintegrable system was viewed as a 
relatively simple and somehow trivial problem JTTJ , but several authors have 
proved that this quantization procedure has to add correction terms to the 
integrals of motion or to the Hamiltonian, these correction terms are of order 
0(h 2 ) [[TJ, The result of the quantum deformation of a superintegrable 



system is realized by the shift of the classical Poisson algebra to some quan- 
tum polynomial associative algebra. The same fact is true in the case of 
quadratic Poisson algebra corresponding to the Yang - Baxter equation ||, 
which is turned to a quantum quadratic associative algebra fLM with four 



generators. The same idea was discussed in reference |10|, where the classi- 
cal problems, which are expressed by a quadratic Poisson algebra are mapped 
to quantum ones described by the corresponding quantum operator quadratic 
algebra. The same shift is indeed true for the superintegrable systems, where 
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the classical ones correspond to the quantum ones and the classical quadratic 
Poisson algebra is mapped to a quadratic associative algebrafT^-[|19|j. 

In this paper we show that the deformation of the classical Poisson algebra 
to a quadratic associative algebra implies a deformation of the parameters 
of the quadratic algebra. The general form of the quadratic algebras, which 
are encountered in the case of the two dimensional quantum superintegrable 
systems, is investigated. In references [0, || [15|, [Tj| [L7], [20|, |21, 



23| was 



conjectured that, the energy eigenvalues correspond to finite dimensional rep- 
resentations of the latent quadratic algebras. Granovskii et al in QTO]] studied 
the representations of the quadratic Askey - Wilson algebras QAW(3). Using 
there the proposed ladder representation, the finite dimensional representa- 
tions are calculated and this method was applied to several superintegrable 

Another method |5], |6| for calculating the fi- 



systems [15, 16, 17, Ell 23 



nite dimensional representations is the use of the deformed oscillator algebra 
[23J and their finite dimensional version which are termed as generalized de- 
formed parafermionic algebras[25|. The main task of this paper is to reduce 



the calculations of eigenvalues to a system of two algebraic equations with 
two parameters to be determined. These equations are universal equations, 
which are valid of all superintegrable systems, with quadratic integrals of 
motion. 

This paper is organized as follows: In section [H] the general form of the 
quadratic Poisson algebra for a two dimensional system with quadratic in- 
tegrals of motion is derived. In section [TT1] the special form of the Poisson 
algebra of the known two dimensional superintegrable systems in flat space is 
written. In section fV| the quantum version of the Poisson quadratic algebra 
is studied. The deformed parafermionic oscillator algebra is reviewed and the 
oscillator realization of the quadratic algebras is realized. The finite dimen- 
sional representations of the quadratic algebras are generated by using the 
technique of deformed parafermionic algebras. The problem is reduced to the 
solution of a system of two algebraic equations in section |VI[ In section [VT i 
the energy eigenvalues of all the known superintegrable systems in the flat 
two dimensional space are determined by solving the appropriate algebraic 
equations. Finally in section V11J there is a discussion of the results of this 
paper. 
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II Quadratic Poisson Algebras 



Let consider a two dimensional superintegrable system. The general form of 
the Hamiltonian is: 

H = a(gi, q 2 )p\ + 2b(q u q 2 )pip 2 + c(q u q 2 )pj + V(q u q 2 ) (1) 

this Hamiltonian is a quadratic form of the momenta. The system is su- 
perintegrable, therefore there are two additional integrals of motion A and 
B. In that section, we consider that, these integrals of motion are quadratic 
functions of the momenta, i.e. they are given by the general forms: 

A= A(q 1 ,q 2 ,p 1 ,p 2 ) =c(q 1 ,q 2 )pl + 2d(q 1 ,q 2 )p 1 p 2 + e(q 1 ,q 2 )pl+ 

+f(Qi, qt)pi + g(qi, Q2)p 2 + Q(qi, q 2 ) 

The integral of motion B is assumed to be indeed a quadratic form, which is 
analogous to above one. 

B= B(q l ,q 2 ,p 1 ,p 2 ) = h(q u l 2 )p\ + 2k(q x , q 2 )pip 2 + l(q u q 2 )p\+ 
+m(q 1 , q 2 )pi + n(q 1 , q 2 )p 2 + S(q 1 , q 2 ) 

By definition the following relations are satisfied: 

{H,A} p = {H,B} p = (2) 

where { . , . } p is the usual Poisson bracket. 

From the integrals of motion A, B, we can construct the integral of mo- 
tion: 

C = {A,B} P (3) 

The integral of motion C is not a new independent integral of motion, which 
is a cubic function of the momenta. The integral C is not independent from 
the integrals H, A and B as it will be shown later. The fact that, the integral 
C is a cubic function of momenta, implies the impossibility of expressing C as 
a polynomial function of the other integrals, which are quadratic functions of 
momenta. Starting from the integral of motion C, we can construct the (non 
independent) integrals {A,C} P and {B,C} p . These integrals are quartic 
functions of the momenta, i.e. functions of fourth order. Therefore these 
integrals could be expressed as quadratic combinations of the integrals H, A, 
and B. Therefore the following relations are assumed to be valid: 

{A, C} p = aA 2 + (3B 2 + 2 7 AB + 5 A + eB + ( (4) 
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and 

{B, C} p = aA 2 + bB 2 + 2cAB + dA + eB + z (5) 

We can take appropriate a linear combination of the integrals A and B and 
we can always consider the case (3 = 0. 

The Jacobi equality for the Poisson brackets induces the relation 

{A,{B,C} p } p = {B,{A,C} p } p 

The following relations 

b = —7, c = —a and e = —5 

must be satisfied. 

The integrals A, B and C satisfy the quadratic Poisson algebra: 

{A,B} p = C 

{A,C} P = aA 2 + 2 7 A, B + 5 A + eB + C (6) 
{B,C} P = aA 2 --fB 2 -2aAB + dA-5B + z 

where a, 7, a are constants and 

S = 5(H) = 5 + 5 X H 

e = e(H) = e + e x H 

( = ((H)=(; + ( 1 H + ( 2 H 2 

d = d(H) = d + diH 

z = z(H) = z + z x E + z 2 H 2 

where 5i, e^, Q, di and Z{ are constants. The associative algebra, whose the 
generators satisfy equations @, is the general form of the closed Poisson 
algebra of the integrals of superintegrable systems with integrals quadratic 
in momenta. 

The quadratic Poisson algebra (|6|) possess a Casimir which is a function 
of momenta of degree 6 and it is given by: 

K= C 2 - 2aA 2 B - 2-/AB 2 - 25 AB- 

-eB 2 - 2(B + laA 3 + dA 2 + 2zA = (7) 
= k + k x E + k 2 H 2 + k 3 H 3 

Obviously 

{K,A} p = {K,B} p = {K,C} p = 
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Therefore the integrals of motion of a superintegrable two dimensional system 
with quadratic integrals of motion close a constrained classical quadratic 
Poisson algebra (|6|), corresponding to a Casimir equal at most to a cubic 
function of the Hamiltonian (|7|). 

In the general case of a superintegrable system the integrals are not nec- 
essarily quadratic functions of the momenta, but rather polynomial functions 
of the momenta. The case of the systems with a quadratic and a cubic inte- 
gral of motion are recently studied by Tsiganov ||26|| . The general form of the 
Poisson algebra of generators A, B, and C is characterized by a polynomial 
function h(A, B), which satisfy the following equations: 

{A,B} p = C 

{A,C} P = dh/OB (8) 
{B,C} P = -dh/dA 

and the Casimir of the algebra is given by 

K = K(H) = C 2 — 2h(A, B), {K, A} p = {K, B} p = (9) 

where h(A, B) is a polynomial function of the integrals of motion A and B. 
In the case of the quadratic Poisson algebra (|5p the form of the function 
h(A, B) is given by equation (|7|): 

h(A,B)= — |v4 3 + aA 2 B + ^AB 2 
-\A 2 + SAB + \B 2 
-zA + CB 

In the general case of a two dimensional superintegrable system, with 
quadratic Hamiltonian, one integral A of order m in momenta and one inte- 
gral B of order n (n > m), the general form of the function h(A, B) can be 
given by the general form: 

h(A, B) = h (A) + h\{A)B + h 2 (A)B 2 

where hi(A) are polynomials of the integrals A and H. The proof of this 
assumption is based on the dependence of the integrals of motion on the 
momenta. For simplicity reasons, the proof of this proposition will not be 
given here. 
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Ill Poisson algebras for superintegrable sys- 
tems 



Let consider the superintegrable systems with quadratic integrals of mo- 
tion, these potentials are given by several authors starting from different 
but comparable points of view. In references |2| the integrals of motion 
are generated by solving the Darboux conditions for integrability of quaratic 
integrals. In [f§]the Hamilton - Jacobi equation is solved by separation of 
variables nad the two dimensional Hamiltonians which are separable in more 
than one coordinate system are obtained. The separation of variables is es- 
sential for solving the quantum counterpart of the superintegrable system 
and the solution of the associate Schrodinger equations is given in ]7| . Using 
this method the quantum superintegrable systems have been solved on the 
sphere^] and the hyperboloid|0, [19| . From classical point of view the super 
integrable are given in ||, while the case of a pseudo Euclidean kinetic term 
has been studied in 0. The extension on the systems with a quadratic and 



a cubic integral of motion is sytematized in |26 . 

In this section we consider the case superintegrable systems gin in ref , 
because in the next sections we study the quantum counterparts of these po- 
tentials. In this paper the following superintegrable systems are considered: 
Potential i): 

2 , 2 , 2 2 , Z 1 ! , ^2 | 
Px +P y + UJ T + ^ + ^J 

This potential has the following independent integrals of motion: 

A = p 2 x + lu 2 x 2 + ^ 
x A 




and 



B = (xp y - yp x f + r 2 (^ + 



The constants, which characterize the corresponding quadratic algebra 
are given by: 

a = -8, 7 = 0, S = 16H, 

e = -16cu 2 , ( = 16(^i + /i 2 V 2 

a = 0, d = 0, z = 16(/i 2 — /-ti V 2 
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the value of the Casimir ([?]) is: 

K = -16 ((// 2 - ^) 2 u 2 + 
Potential ii): 

This potential has the following independent integrals of motion: 

A = pi + 4uj 2 x 2 

and 

£ = {xpy - yp x ) p y + — -ujxy 

The constants, which characterize the corresponding quadratic algebra 
are given by: 

a = 0, 7 = 0, 5 = 0, 

e = -16cj 2 , C = 0, 

a = -6, d = 16if, z = 8/i 2 cu 2 - 8# 2 
the value of the Casimir (0) is: 

K = 

Potential iii): 

2 \ y r r \r + x r — x / y 
This potential has the following independent integrals of motion: 

A = (xp y - yp x ) 2 + r ( + ) 

\r + x r — xj 

and 

H\r — x /i 2 r + x fcx 

£ = xp y - - 7 : — + 7, + TT 

2r r + x 2r r — x 2r 
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The constants, which characterize the corresponding quadratic algebra 
are given by: 

a = 0, 7 = -2, 5 = 0, 

e = 0, C = k(jii - /i 2 ), 

a = 0, d=-8H, z = 4(/ix + /i 2 )# - k 2 /2 

the value of the Casimir (0) is: 

X = 20i - ^a) 2 ^ - A; 2 (/ii + /i 2 ) 

Potential iv): 

1 / 2 n & V r + ^ V r _ X \ 
H = ~ \Px+Py + - + +/^2 

2 \ y r r r J 

This potential has the following independent integrals of motion: 

liAr — x)\/r + u /i 2 (r + x)y/r — u kx 
A = (y Px - xp y ) Py + 

and 

, \ HiXy/r - u il 2 x\Jt + u ky 
B = (x Py - yp x )p x + — 

The constants, which characterize the corresponding quadratic algebra 
are given by: 

a = 0, 7 = 0, 5 = 0, 

e = 2H, C = -/*W2, 

a = 0, d= -2H, z = 

the value of the Casimir (|7|) is: 

K = -k 2 H/2 - k(/4 + /i 2 )/4 



IV The quadratic associative algebra 

The quantum counterparts of the classical systems, which have been studied 
in section |TT|, are quantum superintegrable systems. The quadratic classical 
Poisson algebra (|6]) possesses a quantum counterpart, which is a quadratic 
associative algebra of operators. The form of the quadratic algebra is similar 
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to the classical Poisson algebra, the involved constants are generally functions 
of h and they should coincide with the classical constants in the case h — > 0. 
Let consider the quadratic associative algebra generated by the generators 
{A, B, C}, which satisfy the commutation relations 

[A,B] = C 

[A, C] = aA 2 + (3B 2 + 1 {A, B} + 5A + eB + ( (10) 
[B,C] = aA 2 + bB 2 + c{A,B} + dA + eB + z 

After rotating the generators A and B, we can always consider the case f3 = 0. 
The Jacobi equality for the commutator induces the relation 

[A, [B,C]} = [B, [A,C]} 

the following relations 

b = —7, c = —a and e = —5 

must be satisfied, and consequently the general form of the quadratic algebra 
([T0| ) can be explicitly written as follows: 



[A,B]=C (11) 

[A,C] = aA 2 + 1 {A,B} + 5A + eB + ( (12) 
[B, C] = aA 2 — "fB 2 — a {A, B} + dA-5B + z (13) 
The Casimir of this algebra is given by: 

K = C 2 -a {A 2 , B} - 7 {A, B 2 } + (try - 5) {A, B} + 

+ ( 7 2 - e)B 2 + ( 7 5 - 2C)5+ (14) 
+ f A 3 + {d + f + a 2 )A 2 + {f + a5 + 2z)A 

another useful form of the Casimir of the algebra is given by: 



K= C 2 + f A 3 -f {A,A,B}~1 {A,B,B} + 
+ (^ + d+ 2 -f)A 2 + {-e + 2£) B 2 + 
+ (S + ^){A,B}+(^ + f + ^ + 2z) A+ 



3 / l ' J 1 V3 1 3 1 3 



(15) 



+ ("f + T- 2 CJ5 ■ 3 3 
where 

{A, B, C} = ABC + ACB + BAC + BCA + CAB + CBA 
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This quadratic algebra has many similarities to the Racah algebra QR{3) , 
which is a special case of the Askey - Wilson algebra QAW(3). The algebra 
(|TT| - [13|) does not coincide with the Racah algebra QR(3), if a 7^ in 



the relation (|13|). Unless this difference between (pTOj) and QR(3) algebra a 
representation theory can be constructed by following the same procedures 



as they were described by Granovskii, Lutzenko and Zhedanov in ref. |[T0| . 
P~5| , [U|. In this paper we shall give a realization of this algebra using the 
deformed oscillator techniques ||24|| . The finite dimensional representations 



of the algebra (|i~0"D will be constructed by constructing a realization of the 
algebra with the generalized parafermionic algebra introduced by QuesnepS]. 



V Deformed Parafermionic Algebra 



Let now consider a realization of the algebra (O - O), by using of the 



deformed oscillator technique, i.e. by using a deformed oscillator algebra |24 
j&^^TVj, which satisfies the 

'AT, b ] ] = b\ [AT, b] = -b, b ] b = $ (AO , bb ] = $ (A" + 1) (16) 

where the function $(x) is a "well behaved" real function which satisfies the 
boundary condition: 



$(0) = 0, and $(x) > for x > 



(17) 



As it is well known p^j this constraint imposes the existence a Fock type 
representation of the deformed oscillator algebra, which is bounded by bellow, 
i.e. there is a Fock basis \n>, n — 0, 1, . . . such that 



Af\n >= n\n > 

tf\n >= [n + l)\n + 1 >, 

6|0 >= 

b\n >= J$~(n)\n - 1 >, n 



n 



0,1,... 



(18) 



1,2,... 



The Fock representation (THf) is bounded by bellow. The generalized de- 



formed algebra given in ref [24] is equivalent to several deformed oscilla- 



tor schemes as the Odaka- Kishi - Kamefuchi unification scheme ||27|| , the 



Beckers- Debergh unification scheme ||28|| , The Fibonacci oscillator [p9|] , for 
a discussion of deformation schemes see |30 
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In the case of nilpotent deformed oscillator algebras, there is a positive 
integer p, such that 

= o, (tf) p+1 = 



the above equations imply that 

$(p+l) = 0, (19) 

In that case the deformed oscillator ( |I6|) has a finite dimensional representa- 
tion, with dimension equal to p+1, this kind of oscillators are called deformed 
parafermion oscillators of order p. 

An interesting property of the deformed parafermionic algebra is that the 
existence of a faithful finite dimensional representation of the algebra implies 
that: 

A" (A" - 1) (A" - 2) ■ • • (A" - p) = (20) 

The above restriction and the constraints flTT|) and (|T9|) imply that the general 
form of the structure function $(7V) has the general form|25|: 



$(A0 = M{p + 1 - A0(a + aiA/" + a 2 M 2 + ■■■ a^M^ 1 ) 

A systematic study and applications of the parafermionic oscillator is given 
in references |25|, [31], |32|, |33| . 

We shall show, that there is a realization of the quadratic algebra , such 
that 

A = A (A/") (21) 
B = b (A") + 6 f p (A") + p (A") b (22) 

where the A[x], b[x] and p(x) are functions, which will be determined. In 
that case (|11|) implies: 

C = [A, B] C = tfAA (A") p{N)-p (A") AA (A") b (23) 

where 

(A") = A (M + 1) -A (A") 

Using equations (PH), ([22]) and ( |T2"D we find: 

[A, C] = [A (A") , tfAA (A") p{N)-p (A") AA (A") 6] = 
= tf(AA(M)fp(M)+p(M)(AA(M)fb = 
= aA 2 + 1 {A,B} + 5A + eB + ( = 

= tf{ 1 {A{N+l)+A{N)) + e)p{N)+ 1 > 

+p (A") ( 7 (A (A" + 1) + A (A")) + e) 6+ 
+aA (A") 2 + 2 7 A (A") 6 (A") + 5 A (A") + eB (A") + C 
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therefore we have the following relations: 



(AA (A/")) 2 = ^ (A(Af + 1) + A (A/")) + e (25) 
aA (AT) 2 + 2 7 A (A/") 6 (A/ - ) + 5 A (Af) + efi (AO + C = (26) 

while the function p (A/ - ) can be arbitrarily determined. In fact this function 
can be fixed, in order to have a polynomial structure function $(sc) for the 
deformed oscillator algebra (|T6|). 

The solutions of equation (p5|) depend on the value of the parameter 7, 
while the function b(Af) is uniquely determined by equation ( p6|) (provided 
that almost one among the parameters 7 or e is not zero). At this stage, the 
cases 7 7^ or 7 = 0, should be treated separately. We can see that: 

Case 1: 7 ^ 

In that case the solutions of equations (|25f) and (^6| ) are given by: 



7 

2 . v ■ -, - ^ 



^(A^) = 7t((A^ + «) 2 -1/4-4 (27) 



q((Ar+ M ) 2 -l/4) ae _5 7 
4 ~+~ 2 7 2 



6 (AT) = 

ae 2 -2(5e7+4 7 2 C JL 



(28) 



4 7 4 ((AT+«)2-l/4) 



Case 2: 7 = 0, e ^ 

The solutions of equations (|25|) and (pB"D are given by: 



A(A0 = \fe{N + u) (29) 



b(Ai ) = -a(N ' + uY- -^{Af ' + «)-- (30) 
The constant u will be determined later. 

Using the above definitions of equations A(Af) and 6(A/"), the left hand 
side and right hand side of equation ( |13"D gives the following equation: 

2 $(A/" + 1) (AA (AT) + l) p{M) - 2 $(A0 (AA (AT - 1) - f ) p(A/" - 1) = 
= a A 2 (A/") - 7& 2 (A0 - 2«A (A/) 6(A0 + dA (A/") - 5&(A0 + 2 

(31) 
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Equation ( |HD gives the following relation: 



K 



= $(JV + 1) ( 7 2 - e - 2 1 A (Af) - AA 2 (A/")) p(Af)+ 

+$(AT) ( 7 2 - e - 2 7 A (AT) - AA 2 (Af - 1)) P (Af - 1)- 
-2aA 2 (AT) 6(A0 + (7 2 - e - 2 7 ^ (A/")) 6 2 (AT)+ 
+2 (cry - 5) A (A/") 6(A/") + ( 7 5 - 2£) &(A0+ 
+ §aA 3 (AT) + (d + |a 7 + a 2 ) A 2 (A/") + 



Equations ( |31| ) and (^) are linear functions of the expressions <E> (AT) and 
$ (Af +1), then the function <3> (AT) can be determined, if the function p(Af) 
is given. The solution of this system, i.e. the function $ (Af) depends on two 
parameters u and K and it is given by the following formulae: 

Case 1: 7 ^ 



$(A0 = -3072 7 6 J fsT(-l + 2(JV + u)) 2 - 

— 48 7 6 (a 2 e — a<5 7 + ae 7 — d 7 2 )- 

•(-3 + 2(JV + u))(-l + 2{M + u))\l + 2(JV + 

+ 7 8 (3a 2 + 4a 7 )(-3 + 2(JV + w)) 2 (-l + 2(JV + u))\l + 2(JV + w)) 2 + 
+768(ae 2 -25e 7 + 4 7 2 C) 2 + 

+32 7 4 (-l + 2{M + u)) 2 {-\ - 12{M + u) + 12{M + u) 2 )- 
•(3a 2 e 2 - 6a<5e 7 + 2ae 2 7 + 25 2 7 2 - 4de 7 2 + 8 7 3 z + 4a 7 2 C)- 
-256 7 2 (-l + 2(A/" + u)) 2 - 

•(3a 2 e 3 - 9a5e 2 7 + ae 3 7 + 65 2 e 7 2 - 3de 2 7 2 + 25 2 7 4 + 
+2de 7 4 + 12e 7 3 ^ - 4 7 5 ^ + 12ae 7 2 C - 125 7 3 C + 4a 7 4 C) 




3 • 2 12 ■ 7 8 (A/" + w)(l + M + u )(l + 2(JV + u)) 2 



and 



(33) 



Case 2: 7 = 0, e + 
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$(A0 = 

_ 1 (_K _ _z S_C i e\ _ 

(3d - - 3a£ + 3 (^) 2 - 6f e + Qa^ - 6^) (Af + u) 
+ \ {a 2 + d - a^~e - 3a^ + (^) 2 + 2a$ \ (Af + uf- 
-\ (3a 2 -ay/e- 3a^) (AT + u) 3 + \a 2 (M + u) A 

(34) 

The above formula is valid for e > 0. 



VI Finite dimensional representations of qua- 
dratic algebras 

Let consider a representation of the quadratic algebra , which is diagonal 
to the generator A and the Casimir K. Using the parafermionic realization 
defined by equations fl2lD and (p2|), we see that this a representation diagonal 



to the parafermionic number operator Af and the Casimir K. The basis of 
a such representation corresponds to the Fock basis of the parafermionic 
oscillator, i.e. the vectors \k, n >, n = 0,1,... of the carrier Fock space 
satisfy the equations 

Af\k, n >= n\k, n >, K\k, n >= k\k, n > 

The structure function (j3~3"D (or respectively (j33|) ) depend on the eigenvalues 
of the of the parafermionic number operator Af and the Casimir K. The 
vectors \k, n > are also eigenvectors of the generator A, i.e. 



A\k, n >= A(k,n)\k, n > 



In the case 7 ^ we find from equation ( B7[ ) 

A{k,n) = 1 -({n + u) 2 -l/A-^ 

In the case 7 = 0, e ^ we find from equation (f29|) 

A(k, n) = \fe {n + u) . 
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If the deformed oscillator corresponds to a deformed Parafermionic oscil- 
lator of order p then the two parameters of the calculation k and u should 
satisfy the constrints ( |iT|) and ( |T9| ) the system: 

$(0,u,k) = 
$(p+l,u,k) = 1 1 



then the parameter u = u(k,p) is a solution of the system of equations (p5|). 

Generally there are many solutions of the above system, but a unitary 
representation of the deformed parafermionic oscillator is restrained by the 
additional restriction 

$(x) > 0, for x—l,2,...,p 

We must point out that the system (|35|) corresponds to a representation with 
dimension equal to p + 1 . 

The proposed method of calculation of the representation of the quadratic 
algebra is an alternative to the method given by Granovskii et al |K], [15]. [L6]. 
[17| and reduces the search of the representations to the solution of a system 
of polynomial equations (|3"5]). Also its is applied to an algebra not included in 
the cases of the algebras, which are treated in the above references. We must 
point out, that there are several papers on the representations of quadratic 
( or generally polynomial algebras) |34], [35], [36], [37], |38|, [59[ f4(J , these algebras 



are deformations of the su(2) or osp(l/2) algebras. The general form of the 
quadratic algebra, which is studied in this paper, is different by definition 
from the deformed versions of su(2) or osp(l/2). 

VII Quadratic algebras for the quantum su- 
perintegrable systems 

In this section, we shall give an example of the calculation of eigenvalues 
of a superintegrable two-dimensional system, by using the methods of the 
previous section. The calculation by an empirical method was performed in 
and the solution of the same problem by using separation of variables was 
studied in . Here in order to show the effects of the quantization procedure 
we don't use H = 1 as it was considered in references and . That means 
that the following commutation relations are taken in consideration: 

[x,p x ]=ifi, [y,p y ] = ih 
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VH-a Potential i) 



, , 1 / 2 2 2 2 /'l /'2 



This potential has the following independent integrals of motion: 



and 



A 2 i 2 2 i /^l 

A = p x + uj x + — 



B = (x Py - yp x f + r 2 f ^ + p 



The constants, which characterize the corresponding quadratic algebra (10). 
are given by: 

a = 8ft 2 , 7 = 0, 5 = -16h 2 H, 

e = i6frV, C = -16^ 2 (yUi + /i 2 )^ 2 + 8^V 

a = 0, = 16^ 4 , z = -lQh 2 (fi 2 - fi^uj 2 - 16fr 4 # 

the value of the Casimir flTj]) is: 

if = 16ft 2 ((^ 2 - /ii) 2 ^ 2 + 4/iiif 2 ) - 16h 4 (3H 2 + 27iV - 2(/i! + // 2 )) 

For simplicity reasons we introduce the positive parameters k% and hi, which 
are related to the potential parameters H\ and fi 2 by the relations: 

A*i = (*? " ^) ^ Ms = (* 2 2 - ^) ^ 
This quadratic algebra corresponds to the case 7 = and e > of the 



algebra given by equations (]TT|— 113[) . In that case, the structure function 
of the deformed parafermionic algebra of Section |V| can be given by the simple 
form: 



$(z)= lQh A (x + u-\-!f) (x + u-\ + ^ 



k2 E_ 

2 1 2 2hui 



In the above formula E is the eigenvalue of the energy. The values of the pa- 
rameters u and E corresponding to the a representation of the parafermionic 



17 



algebra of dimension p + 1 are determined by the restrictions (^), which are 
transcribed as: 

$(0) = 0, $(p+l) = 

One should notice, that only the solutions which correspond to positive eigen- 
values of the integral A must be retained. The acceptable solutions are four 
and correspond to the following values of the parameters u and E: 

U= 2 + — E = ^{P + 1 + 2 ) 

where = ±1. The corresponding structure function is 

<3>(x) = 16h 4 x (p + 1 — x) (x + eiki) (p + 1 — x + e-ik-i) 
The corresponding eigenvalues of the operator A are given by: 

Aim) = Ahu H — 1 1 ^ 2 ^ ra = 0,1, ... ,p 

The structure function should be a positive function, for x — 1, 2, . . . ,p 
therefore the constants fc« are restricted by the relations: 

eih > -1, e 2 k 2 > -1 

VH-b Potential ii) 



H = \ [vl + P 2 y + (Ax 2 + 



,2 



This potential has the following independent integrals of motion: 

A = pi + Auj 2 x 2 



and 

B = - {xp y - yp x , p y \ + — - u xy 

^ y 



The constants, which characterize the corresponding quadratic algebra fllPf), 
are given by: 

a = 0, 7 = 0, 5 = 0, 

e = 16/iV, C = 0, 

a = 6h 2 , d=-16h 2 H, z = -8h 2 (fxco 2 -H' 2 ) + 6h 4 uj 2 
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the value of the Casimir ( fH| ) is: 

K = Q4h 4 u 2 H 

For simplicity reasons we introduce the positive parameters k, which is re- 
lated to the potential parameter fi by the relation: 

„ = (*" - i) h 2 

This quadratic algebra corresponds to the case 7 = and e > of the 



algebra given by equations ( ]TTH l3l ) - In that case, the structure function 
of the deformed parafermionic algebra of Section |V| can be given by the simple 
form: 

, / 1\ / 1 k E \ ( I k E \ 

•(*) = !»*<■< (* + «-jJ ^ + "-2-2-2&:j \ X + n -2 + 2-^) 

In the above formula E is the eigenvalue of the energy. The values of the pa- 
rameters u and E corresponding to the a representation of the parafermionic 
algebra of dimension p + 1 are determined by the restrictions (0), which are 
transcribed as: 

$(0) = 0, $o+i) = o 

One should notice, that only the solutions which correspond to positive eigen- 
values of the integral A must be retained. The acceptable solutions are four 
and correspond to the following values of the parameters u and E: 

1 ( ek\ 

u=-, £ = 2fi^ P+ l + -j 

where e = ±1. The corresponding structure function is 

= Ah 3 x (p+l—x)(p + l—x + ek) 

The structure function should be a positive function, therefore the values of 
the parameter k are restrained by 

ek > -1 

The eigenvalues of the operator A are given by: 

A(m) = 4hu)(m + -), m = 0, 1, . . . ,p 
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VII-c Potential iii) 

2 \ y r r \r + x r — xj J 
In ref |7j the parabolic coordinates have been used: 

y = h, p y = ef^Pt + ewP^ 

[£,Pt]=ih, [r],p v ] = ih 

For comparison reasons we quote all the relations in both, cartesian and 
parabolic systems, so 

This potential has the following independent integrals of motion: 
and 

B = \ {{xPy - yP^Py} + 

The constants, which characterize the corresponding quadratic algebra flIU|) , 
are given by: 

a = 0, 7 = 2h 2 , 5 = 0, 

e = -ft 4 , ( = -h 2 k([ii - fi 2 ), 

a = 0, d = 8h 2 H, z = -h 2 (4(/i! + fi 2 )H - k 2 /2) + Ti A H 

the value of the Casimir (|14|) is: 

K = -h 2 (2(/ix - fi 2 ) 2 H - k 2 (fn + fMi)) - 2^ 4 (V + fMi)H n&H 

For simplicity reasons we introduce the positive parameters k% and ki, which 
are related to the potential parameters fix and fi 2 by the relations: 




This quadratic algebra corresponds to the case 7 7^ of the algebra 
given by equations (pHllf). In that case, the structure function (|33"D of the 
deformed parafermionic algebra of Section [V] can be given by the simple form: 

$0) = 3 ■ 2 14 fr 16 (2x-l + k 1 + k 2 ) (2x - 1 + ifej - k 2 ) {2x - 1 - ifex + fc 2 ) • 
■ (2x - 1 - ifei - fc 2 ) (8k 2 Hx 2 - 8h 2 Hx + 2h 2 H + k 2 ) 

In the above formula E is the eigenvalue of the energy. The values of the pa- 
rameters u and E corresponding to the a representation of the parafermionic 
algebra of dimension p + 1 are determined by the restrictions ([^), which are 
transcribed as: 

$(0) = 0, $O+l) = 

One should notice, that only the solutions which correspond to positive eigen- 
values of the integral A must be retained. The acceptable solutions are four 
and correspond to the following values of the parameters u and E: 

1 k 2 

u — — (2 + e\k\ + €2^2) , E = ^ 77 

2 l >' 2h 2 (2(p+l) + e 1 k 1 + e 2 k 2 ) 2 

where = ±1. The corresponding structure function is 

= 3 - 2 20 -k 2 h w -x{p + 1 - x) (x + eiifei) (x + e 2 k 2 ) ■ 
(x + £l kx + e 2 k 2 ) {2{p+1)+eiki+e2k2f 

The eigenvalues of the operator A are given by the formula: 

A(m) = h 2 [m + eiki + e 2 k 2 + - J , m = 0, 1, . . . , p 

The positive sign of the structure function for x = l,2,...,p is obtained 
when: 

eiki > —1, £2^2 > — 1, and e\ki + e 2 k 2 > —1 
VH-d Potential iv) 



v' — x 



= &w (I (pI + Py)+ k + M + w 
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This potential has the following independent integrals of motion: 

A= \ ({(y Px - xp y ),py} + aO^±H - Mr±^H _ M) = 
= - ^ + fc fa 2 - £ 2 ) + 2 ^ - M)) 

and 

The constants, which characterize the corresponding quadratic algebra (|H)D, 
are given by: 

a = 0, 7 = 0, 5 = 0, 

e = -2ft 2 #, C = fr 2 /Ui/i 2 /2, 

a = 0, d = 2h 2 H, z = -h 2 (fij - nl) /4 

the value of the Casimir (|14|) is: 

= h 2 k 2 H/2 + fr 2 A;(/i 2 + fi 2 2 ) /4 + fr 4 # 2 

This quadratic algebra corresponds to the case 7 = and e > of the 
algebra given by equations (11 -|l3l). It is worth noticing that the algebra is 
extremely simple, which can be reduced to the usual su(2). We prefer to treat 
this algebra with the proposed methods in this paper for pedagogical reasons. 
The existence of the finite dimensional representations of this algebra implies 
that, the coefficient e in equation flTjD should be positive, therefore the energy 
operator H must have energy eigenvalues E < 0. For simplicity reasons we 
introduce the new parameters: 

e = y/=2E/h, X = k/h 2 , 

u 1 = fi 1 /h 2 , v 2 = n 2 /h 2 , v 2 = v\ + v\ 

The structure function ( [34"1) of the deformed parafermionic algebra of 
Section [V] can be given by the form: 

® {x) = i£ " X£2 +2 ( x+u ~ \> 3 ) - Xe2 - 2(* + « - \y 

In the above formula the parameter e is related to the the eigenvalue E of 
the energy. The values of the parameters u and e, corresponding to the a 
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representation of the parafermionic algebra of dimension p+1, are determined 
by the restrictions (|35|), which are transcribed as: 



$(0) = 0, $(p+l) = 

The first condition can be used for determining the acceptable values of the 
parameter u. Two possible solutions are found: 

v\ - Xe 2 + e 3 
u = u x = — 3 (36) 

v\ - Xe 2 - e 3 

u = u 2 = — z (37) 

Using these solutions and the condition Phi(p + 1) = 0, we find that the e 
must satisfy two possible cubic equations: 

Ui — ► 2(p + l)e 3 - 2Xe 2 + v 2 = (38) 
u 2 — ► 2(p + l)e 3 + 2Xe 2 - v 2 = (39) 



If £ is a solution of equation (|38| ) then — e is the solution of the other equation 
(l3"9f) , therefore there is almost one solution which is positive. This solution 
leads to the structure function: 



s 2 



$(x) = —x (p + 1 — x) 



which is positive for x = 1, 2, ... ,p. 



VIII Discussion 

If we compare the quadratic associative algebra, introduced in section [TV] 
with the corresponding Poisson algebra given in section |J we see that in 
general the quantum constants are similar to the classical ones up to an factor 
equal to — h 2 , but there are quantum corrections of order h A and h 6 . The 
knowledge of the classical constants of the Poisson algebra is not sufficient to 
reproduce the rules of quantum associative operator algebra. Therefore, the 
passage from the classic Poisson algebra to the non commutative quantum 
algebra can not be realized by simple replacements of the Poisson brackets 
by commutators and by a symmetrization procedure. 
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The energy eigenvalues of section |V11| corroborate the results of reference 
|7j (the differences in the case of the potential iv are due to some misprints 
in that reference). The calculation of the energy eigenvalues in reference 
|7j was achieved by solving the corresponding Schroedinger differential equa- 
tions, while in this paper the energy eigenvalues are obtained by algebraic 
methods. The advantage of the proposed method is that, the energy eigen- 
values are reduced to simple algebraic calculations of the roots of polynomial 
equations, whose the form is universally determined by the the structure 
functions (]33f) , (|33|) and the system (]35|) . These equations are valid for any 
two dimensional superintegrable system with integrals of motion, which are 
quadratic functions of the momenta. The same equations should be valid 
in the case of two dimensional superintegrable systems in curved space HT 



The superintegrable systems bring up for discussion the open problem of the 
quantization of a Poisson algebra in a well determined context, because these 
systems and their quantum counterparts are explicitly known. 
From the above discussion several open problems are risen: 

• The calculation of the classical Poisson algebras and their quantum 
counterparts for the totality of the two dimensional problems in curved 
space. This study will lead to the calculation of the energy eigenvalues 
by algebraic methods. 

• The two dimensional superintegrable systems are classified by the val- 
ues of the constants of the Darboux conditions ||. The relation of 
these constants with the constants of the quadratic Poisson algebra is 
not yet known. 

• The Poisson algebras for the Drach superintegrable systems with a 
cubic integral of motion were written by using a classical analogue of 
the deformed parafermionic algebra p6| . Their quantum counterparts 



and the calculation of their energy eigenvalues are is a topic under 
investigation. 

The Poisson algebras and the associated quantum counterparts for the 
three dimensional superintegrable systems are not yet fully studied. 



Recently |42| the quantum quadratic algebras have been written down, 
but a systematic calculation of energy eigenvalues was not yet per- 
formed. 
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The above points show that, the study of non linear Poisson algebras and 
their quantum counterparts is a topic of interest. 
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